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Abstract 

In this paper we study the optimal dividend problem for a company whose surplus 
process evolves as a spectrally positive Levy process. This model including the dual model 
of the classical risk model and the dual model with diffusion as special cases. We assume 
that dividends are paid to the shareholders according to admissible strategy whose divi- 
dend rate is bounded by a constant. The objective is to find a dividend policy so as to 
maximize the expected discounted value of dividends which are paid to the shareholders 
until the company is ruined. We show that the optimal dividend strategy is formed by a 
threshold strategy. 
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1 INTRODUCTION 



Recently, dividend optimization problems for financial and insurance corporations have 
attracted extensive attention. How should corporation pay dividends to its shareholders? 
A possible goal is that the company tries to maximize the expectation of the discounted 
dividends until possible ruin of the company. In recent years, quite a few interesting papers 
deal with the optimal dividend problem in the dual model. For example, Avanzi et al. 
(2007) consider the model which is dual to the classical risk mode, the authors studied 
how the expectation of the discounted dividends until ruin can be calculated when gain 
distribution has a exponential distribution or mixtures of exponential distributions and 
shows how the exact value of the optimal dividend barrier can be determined, if the jump 
distribution is known. Avanzi and Gerber (2008) examined the same problem for the dual 
model perturbed by diffusion. Moreover, they pointed out that "the optimal dividend 
strategy in the dual model is a barrier strategy. A direct proof that optimal strategy is a 
barrier strategy is of some interest but has not been given to our knowledge; the proof in 
Bayraktar and Egami (2008) is for exponential gains only." Yao et al. (2010) considered 
the optimal problem with dividend payments and issuance of equity in a dual risk model 
without a diffusion, assuming proportional transaction costs, they found optimal strategy 
which maximizes the expected present value of the dividends payout minus the discounted 
costs of issuing new equity before ruin. In addition, for exponentially distributed jump 
sizes, closed form solutions are obtained. Dai, et al. (2010, 2011) considered the same 
problem as in Yao et al. (2010) for a dual risk model with a diffusion with bounded gains 
and exponential gains, respectively. Avanzi, Shen and Wong (2011) determined an explicit 
form for the value function in the dual model with diffusion when the gains distribution is a 
mixture of exponentials. They showed that a barrier dividend strategy is also optimal and 
conjectured that the optimal dividend strategy in the dual model with diffusion should 
be the barrier strategy, regardless of the gains distribution. Until recently, Bayraktar, 
Kyprianou and Yamazaki (2013) using the fluctuation theory of spectrally positive Levy 
processes, show the optimality of barrier strategies for all such Levy processes. 
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Barrier strategies often serve as candidates for the optimal strategy when the dividend 
rate is unrestricted. However, if a barrier strategy is applied, ultimate ruin of the company 
is certain. In many circumstances this is not desirable. This consideration leads us 
to impose restriction on the dividend stream. Ng (2009) considered the dual of the 
compound Poisson model under a threshold dividend strategy and derived a set of two 
integro-differential equations satisfied by the expected total discounted dividends until 
ruin and showed how the equations can be solved by using only one of the two integro- 
differential equations. The cases where profits follow an exponential or a mixture of 
exponential distributions are then solved and the discussion for the case of a general 
profit distribution follows by the use of Laplace transforms. He illustrated how the optimal 
threshold level that maximizes the expected total discounted dividends until ruin can be 
obtained. In this paper we provide a uniform mathematical framework to analyze the 
optimal control problem with dividends payout for a general spectrally positive Levy 
process. We assume that dividends are paid to the shareholders according to admissible 
strategies whose dividend rate is bounded by a constant. Under this additional constraint, 
we show that the optimal dividend strategy is formed by a threshold strategy. 

The rest of the paper is organized as follows. Section 2 presents the model and for- 
mulates the dividend optimization problem. Section 3 discusses the threshold strategies. 
Explicit expressions for the expected discounted value of dividend payments are obtained, 
and in Section 4 we give the main results, it is shown that the optimal dividend strategy 
is formed by a threshold strategy with parameters b* and a under the constraint that 
only dividend strategies with dividend rate bounded by a are admissible. 

2 The model and the optimization problem 

Let X = {X t }t>o be a spectrally positive Levy process with non-monotone paths on a 
filtered probability space (ft, J 7 , F, P), where F = (J r t)t>o is generated by the process X 
and satisfies the usual conditions. The Levy triplet of X is given by (c, a, n), where c > 0, 
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a > and II is a measure on (0, oo) satisfying 

(1 Ax 2 )U(dx) < oo. 



o 



Denote by P x for the law of X when X = x. Let E x be the expectation associated with 
P x . For short, we write P and E when X = 0. The Laplace exponent of X is given by 

1 1 f°° 

*(9) = -logEe- eXt = c9 + -a 2 9 2 + / ( e - fe - 1 + 9xl {]x]<1} )U(dx), (2.1) 
£ * Jo 

where 1a is the indicator function of a set A. In the sequel, we assume that — ^'(0+) = 

E(Xl) > which implies the process X drifts to +oo. It is well known that if yU(dy) < 

oo, then E(X X ) < oo, and E(Xl) = — c + yll(dy). Note that X has paths of bounded 

variation if and only if 



POO 

(7 = and / (1 A x)U(dx) < oo. 
Jo 



In this case, the Laplace exponent (2.1) can be written as 



/•oo 

1 r 2n2 | / 



= c 9 + -a 2 ^ + / (e- 0x - l)n(tfa), (2.2) 



o 



with Co = c + Jg 1 xll(da;) the so-called drift of X. 

For an arbitrary spectrally positive Levy process, the Laplace exponent ^ is strictly 
convex and lim^oo ^(9) = oo. Thus there exists a function $ : [0, oo) — > [0, oo) defined 
by $(g) = sup{# > O\V(0) = q} such that = g, g > 0. 

For more details on spectrally positive Levy processes, the reader is referred to Bertoin 
(1996) and Kyprianou (2006). 

Assume the canonical decomposition of X is given by 

X t = -ct + oB t + J t ,t> 0, (2.3) 

where {B t ,t > 0} is a standard Wiener process, {Jt,t > 0} is a pure upward jump 
Levy process that is independent of {B t ,t > 0}. In addition J = 0. Note that the 
dual model with diffusion in Avanzi and Gerber (2008) corresponds to the case in which 
U(dx) = XF(dx), where A > is the Poisson parameter and F is the distribution of 



individual gains, and the rate of expenses is given by Cq = c + Jq xll(dx). In particular, 
when a — 0, the model deduces the so-called dual model in Avanzi, Gerber and Shiu 
(2007). 

We now recall the definition of the q— scale function W^ q \ For each q > there exits a 
continuous and increasing function : R — > [0, oo), called the g-scale function defined 
in such a way that W^ q \x) = for all x < and on [0, oo) its Laplace transform is given 
by 

poo 1 

/ e- 8x W^(x)dx = — , > <%). (2.4) 

Jo w\p) - q 

We will frequently use the following functions 



Z {q \x) = l + q f W (q) (y)dy, 
Jo 



x G 



and 



Z {q \x) = f* Z( q \z)dz, x e 
Jo 



Note that 

Z iq \x) = 1, Z iq \x) =x, x<0. 

The following facts about the scale functions are taken from Chan, Kyprianou and Savov 
(2011). If X has paths of bounded variation then, for all q > 0, W^|(o,oo) £ C^O, oo) 
if and only if II has no atoms. In the case that X has paths of unbounded variation, it 
is known that, for all q > 0, VF (9) |( 0jOO ) G C^O, oo). Moreover if a > then C^^oo) 
may be replaced by C 2 (0, oo). Further, if the Levy measure has a density, then the scale 
functions are always differentiable. It is well known that 



and 



— , if X has paths of bounded variation, 



CO 



0, otherwise, 



4, if <r ± 0, 



W i5) '(0+) = { q+nff- 00 ) , if X is compound Poisson 

1 Cq 

oo, if a = and 11(0, oo) = oo. 

In all cases, if E{X^ > 0, then W^(oo) = l/E{X 1 ). Ifq > 0, then W®(x) ~ e^ x /^\^(q)) 
as x — > oo. 
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We assume that the surplus process of the company is modeled by (2.3) if no dividends 
are paid. An admissible (dividend) strategy n = {LJ\t > 0} is given by a nondecreasing, 
right-continuous and F-adapted process starting at 0. Let U 7 " = {U£ : t > 0} be the 
company's surplus, net of dividend payments, at time t. Thus, 



In this article we are interested in the case that tt only admits absolutely continuous 
strategies such that dDJ = l^fydt, l w (t) satisfies < l n (t) < a, where a is a ceiling 
rate. Moreover, we make the assumption that, if X has paths of bounded variation, 



where q > is an interest force for the calculation of the present value and T n is the time 
of ruin which is defined by 



We denote by S the set of all the admissible dividend strategies. The objective is to solve 
the following stochastic control problem: the maximal dividend-value function, which is 
defined as 



and to find an optimal policy 7r* e S that satisfies V(x) = V n *(x) for all x > 0. In this 
paper, we will prove that the optimal dividend strategy is formed by a threshold strategy 
with parameters b* (the definition of b* is given by (4.4)) and a: whenever the controlled 
risk process is below b*, no dividends are paid; however, when the controlled risk process 
is above this level, dividends are paid continuously at the maximal admissible rate a. 



U? = X(t)-L«, t>0. 



a < c + J Q yTl(dy). We define the dividend- value function V n by 




r n = inf{t > 0\U? = 0}. 




(2.5) 



3 Threshold dividend strategies 



In this section, we assume that the company pays dividends according to the following 
threshold strategy governed by parameters b > and a > 0. Whenever the modified 
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surplus is below the threshold level b, no dividends are paid. However, when the surplus 
is above this threshold level, dividends are paid continuously at a constant rate a that does 
not exceed the rate of expense c. We define the modified risk process U b = {Ub(t) : t > 0} 
by U b (t) =X t - D b (t), where D b (t) = a f Q 1 {u b (t)>b}dt. Let D b denote the present value 
of all dividends until time of ruin T, 

D b = a e~ qt l {Ub{t)>b} dt 
Jo 

where T = inf{t > : U b (t) = 0} with T = oo if U b (t) > for all t > 0. Here q > is the 
discount factor. Denote by V(x, b) the expected discounted value of dividend payments, 
that is, 

V(x,b) = E(D b \U b (0) = x). 
Clearly, < V(x,b) < | and lim,,^ V(x, b) = |. 

Denote by T the extended generator of the process X, which acts on C 2 function g 
defined by 

1 f°° 
Ag(x) = -a 2 g"{x) - cg'(x) + J [g(x + y) - g(x) - g'(x)yl {lyl<1} ]U(dy). (3.1) 

Theorem 3.1. Assume that, as a function of x, V(x, b) is bounded and twice continuously 
differentiable on (0, b) U {b, oo) with a bounded first derivative. Then V(x, b) satisfies the 
following integro- differential equations 

AV(x, b) = qV(x, b), 0<x<b, (3.2) 

AV(x, b) - aV'(x, b) = qV(x, b) - a, x > b, (3.3) 
and initial condition V(0,b) = 0. 

Proof Applying Ito's formula for semimartingales one has 
E x [e-^V(U b (tAT),b)] =V(x,b) 

+E X / e-* 8 [(A-q)V(U b (s),b) - al {Ub{s)>b} V'(U b (s),b)}ds. 
Jo 

Letting t — >■ oo and note that V(0, b) = we have 

V(x, b) = aE x [ e- qt l {Ub(t)>b} dt 
Jo 
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if and only if 

(A - q)V(x, b) - al {x>b} V'(x, b) = -al {x>b} . 

This ends the proof. 

Remark 3.1. From (3.2) and (3.3) one can prove that V satisfying the following 

V(b-,b) = V(b+,b)=V(b,b), (3.4) 
(ci + a)V'(b+, b) - dV'(b-, b) = a, if a = 0, (3.5) 
V'(b+, b) = V'(b-, b), if a ^ 0. (3.6) 

Define the first passage times, with the convention inf = oo, 

T+ = inf{i > : U b (t) > b}, T~ = inf{t > : U b (t) < b}. 

For q > 0, let 

$i(g) = sup{# > 0\V(6) + a9 = q}. 

The following result generalized the result of Ng (2009, Theorem 2) in which only the 
dual of the classical insurance risk model was considered. 

Theorem 3.2. For x > b, 

V(x, b) = ^+ (v(b, b) - ^ e -*i(«)(*-*). (3.7) 

Proof By using the strong Markov property of U b at T b as in Yin et al. (2013), we 
have 



V(x,b) = *-^E x (e-^,T-<oo) 

+E x (e-^~V(U b (T b -),b),T b - < oo). 

The result (3.4) follows since P x (U b (T b ) = b) = 1 and 

E x (e- qT > ,T~ < oo) = exp(-$i(g)(a; - b)). 
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Theorem 3.3. For < x < b, 

V(x,b) = -\a 2 W {q) {b-x) 



+ V(b, b) ^ e *i(«)( 6 -*) - ia 2 $ 1 (g)^W(6 - x) 

+a<5> 1 (q)e^ (b - x ^ f * W {q \z)e~^ q)z dz 
Jo ; 

+ ^ (^a 2 ^ 1 (q)W^(b -x)- e^^^iq) jf * Z( q \z)e~^ z dz 

+ -e* 1 ^^) (g -a$i(g)) T "V^e"* 1 * 9 ^*. (3.8) 
Q Jo 



where 



Here 



\a 2 W^{b)-B{b)-C{h) 

V[h ' h) = Mb) 



A(b) = e* 1 *^ - jU 2 $i(g)W^(&) + a$ 1 (g)e* l(,?)(6) / W^>(z) e -* l(9) *efc, 
2 Jo 

B(b) = ^ Q(7 2 $i(g)W (9) (&) - e $l(?)(k) $i(?) Z (q \z)e~^ q)z dz 

pb 

C(b) = -e* l(?)(6) (g - a$i(g)) / W (9) (z) e -* l(9) *dz. 
9 Jo 



In particular, when a = 0, 

pb — x 



V(x,b) = V(b, 6)e* 1 ^( 6 - x ^l + a$ 1 (q) W {q \z)e~^ q)z dz 
+ - e ^)(b-x) (g _ ^ X W (q \z)e-^ z dz 

Q JO 

- -e* 1 ^ 6 -^^) f X Z^(z)e-^ z dz, (3.9) 
9 Jo 



where 

V(b,b) 



^i(g) J* Z (q \z)e-^ (q >dz - f (g - a$i(g)) {fl^je-* 1 ^ 2 ^ 
1 + a$i(g) / 6 W(9)(z)e-*i(9) z dz 



Proof We first assume that A := J* °°n(dy) < oo. Set c = c + J* yU(dy) and 
f(y)dy = then / is a probability density on (0, oo). In this case, in view of (3.2), 
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the integro-differential equation (3.7) can be written as 
1 



a 2 V"(x, b) - c V'(x, b) 



-X 



V(b, b) - 



a 



Q 



exp(-$i(g)(x + y- b))f(y)dy 



b—x 



-X [ V(x + y,b)f(y)dy-—(l-F(b-.,)) 
Jo 1 

+ (A + q)V(x,b), < x < b, 



Xa 



(3.10) 



where F is the distribution function of /. Replace the variable x by z = b — x, and define 
W by W(z, b) = V(b - z, b), < z < b. The (3.10) becomes 

1 



a 2 W"(z, b) + c W'(z, b) 



-A 



W(0,b) - 



a 



exp(-$i(g)(y - z))f{y)dy 



Xa 



-A / W(y,b))f(z-y)dy (1 - F(z)) 

Jo Q 

+ (\ + q)W(z,b), 0<z<b, 



(3.11) 



with initial condition W(0, b) = V(b, b) and boundary condition W(b,b) = 0. We extend 
the definition of W by (3.11) to z > and denote the resulting function by w. Then we 
have 



-a 2 w"(z) + c w'(z) = -A 



w(0) - - 
1 



exp(-$i(g)(y - z))f{y)dy 
Xa 



-X f Z W (y))f( Z -y)dy-—(l-F(z)) 

Jo Q 

+(X + q)w(z), z>0, 



(3.12) 



with w(0) = V(b, b) and w(b) = 0. 



For a function g, denoted by g the Laplace transform of g, i.e. g(£) = / °° e ^ y g{y)dy. 
Then the Laplace transform w for w can be easily determined from Eq. (3.11) as 

\a\-\ + Zw(0)) + c o w(0) + - 1) - p^y^W - f )(/(<&i(?)) - /(0 



(0 



1^2 + ^ - (A + g) + A/(0 



(3.13) 



Note that 



POO 

/ e^W {q \x)dx = 

Jo t;V 2 £ 2 



e + c Z-(X + q) + Xf(0 



poo px 

/ e-*dx / W<*\y)dy= — 

Jo Jo £(o°" 



1 



((V( 2 + c„e-(A + ?) + A/(())' 



e(ve + c £-(A + g) + A/(o) 



/ e-^dx dy W {q \z)dz 
Jo Jo Jo 

Now inverting (3.13) gives 

w(z) = - 1 -a 2 W^(z)+ 1 -a 2 w(0) £ W®(z - y)6' (y)dy 

+c w(0)W^(z) + \w(0) (\V iq) * (/ - /($ i(g))5 ) * l) (z) 

\a ( Z (q) - 1 \ 
+y ^— *(*> -/)**) (z) 

(3-14) 

where 5 is the Dirac delta function at 0, hi * h 2 stands for convolution of hi and h 2 and 
l(z) = exp(<&i(q)z). After some tedious calculations, we get 

[ Z W^(z-y)5' (y)dy = W^(z), 
Jo 

A (w {q) * (/ - f($i(q))6 ) * ij (z) = a$i(q) J' W^(z -y)e^dy 

+e *i(5)* _ I^W^'^) 

-(l* 2 ®i(q) + co)wM(z), 

\(ZM-l)*(6 -f)*l)(z) = (^<T 2 q$i(q) + c q^ £ W^(z - y)e*^dy 

+ \a 2 qW {q \z) -q f Z iq \z - y)e^ y dy. 

Substituting the three expressions above into (3.14) we arrive at 

w(z) = -\a 2 W^\z) 

+ w(0) ^e*i(«> - ^^(^(''(z) 

+ ^ Qfj^^W^ - e* 1 ^^) jT Z^(y)e-^ y dy^j 

+ -e*^ z (q - a$i(q)) / W^e-* 1 ^^, (3.15) 
Q Jo 

and the result (3.8) follows since V(x, b) = w(b — x) and w(0) = V(b, b). 
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Now, we assume that A := / °° H(dy) = oo. Let IT n be measures on (0, oo) defined by 

U n (dx) = n(dx)l {( i j0o)} (x), n > 1. 

Then we have 

n n (cfo) < n 2 / x 2 n(rfa;) + / (1 A x 2 )U(dx) < oo. 
Set c„ = c + yll n (dy) and f n (y)dy = ; then for each n > 1, /„ is a probability 



density on (0, oo). Similar to (3.10) we consider the following integro-differential equation 

exp(-O n (g)(a; + y- b))f n {y)dy 



^a 2 V n "(x,b)-c n V n '(x,b) = -A 



K(M)-- 

6 * V n (x + y, b)f n {y)dy - ^ (1 - F n (6 - x)) 



-v 

Jo 

+ (A + q)V n (x, b), 0<x <b, (3.16) 



where Q„(q) = sup{6> > 0\V n (d) + a9 = q}. Here, 

1 f°° 

*n(0) = ^ + -a 2 ^ 2 + y (e- fe - 1 + 9xl {lxl<1} )U n (dx). (3.17) 
Repeating the same argument as the case that A := f£°IL(dy) < oo, we obtain 
V n (x,b) = - l -o 2 W^\b-x) 

+ V n (b, b) ^ e e-W(»-*) - i (T 2 n (g)H/W(6 - x) 

+ae n (g)e " (9)( ^ :E) ^ ' ' W^\z)e- 6n( - q)z dz^j 

+ ^ (^a 2 Q n (q)W^(b -x)- e e ^ b -^Q n (q) J* * \z)e~ e ^ z dz^j 

+ - e e »MM(,_ a e n (g)) T Vi^zje- 8 "^, (3.18) 
9 Jo 

where H^i 9 ' 1 and are scale functions corresponding to the process X n with Laplace 
exponent \l/ n . Since lim^oo ^ n (9) = ^(9), then X n — > X weakly as n — > oo. Thus 
lim^ooK = V, lim^ooWi^ = W {q \ Xim^^Z^ = Z (<?) and lim,,^ 6„ = $i. Conse- 
quently, (3.8) still holds for this case. This ends the proof of Theorem 3.3. 

Remark 3.2. For a > 0, from the graph of 

+ a$i(g) = q 
12 



one can verify that $i(g) — > when a — > oo. After some simple calculations we get 
limo-foo ot$i(q) = q, 

lim a(q - a$i(<?)) = qV'(0+), 

a^oo 

and 

lim V(b 6) = I> + 1^±)_^±) 

As a result, for < x < b, we arrive at 

M 



which is the expected discounted value of dividend payments for the barrier strategy. See 
Lemma 2.1 of Bayraktar, Kyprianou and Yamazaki (2013). 

4 Optimal dividend strategy 

Suppose the maximal dividend- value function V is C 2 (0, oo) (resp. C 1 (0, oo) when X is 
of unbounded (resp. bounded) variation. Standard Markovian arguments, see Fleming 
and Soner(1993), formally yield that V(x) satisfies the following Hamilton- Jacobi-Bellman 
(HJB) equation 

TV(x) - qV{x) + sup {r(l - V\x))} = (4.1) 

0<r<a 

with V(0) = 0, where T is the extended generator of the process X, which defined by 

1 f°° 
Tg(x) = -a 2 g"{x) - cg'(x) + J [g(x + y) - g(x) - g'(x)yl {lyl<1} ]U(dy). (4.2) 

It follows from the HJB equation (4.1) that an optimal dividend policy has to fulfil r = 
if V'(Uf_) > 1; r = a if V'(U£_) < 1. In some situations the optimal dividend strategy 
is a threshold strategy. In fact if V'(x, 0) < 1 for x > 0, then the threshold strategy 
with b* = is optimal, if V'(x, b*) > 1 for x < b* and V'(x, b*) < 1 for x > b*, then the 
threshold strategy with b* > is optimal. The optimal threshold b* can be obtained by 
V'(b*, b*) = 1. From those facts one sees that if V{x, b*) is a concave function on (0, oo), 
then the optimal dividend strategy is a threshold strategy. 
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It follows from (3.7) that 

V'(x, b) = -$i(g) (v(b, b) - ^ e ~^)i^) i x> b. (4.3) 

If $i(g)f < 1, then V'(x,0) < 1 since ^(0,0) = 0. Thus b* = 0. Now, suppose that 
$i(g)| > 1, we get the condition for b*\ 

-Sxte) (v(b*,b*)-^}=l. 
Or, equivalently, b* is the solution of the equation 



From (3.7) and (3.8) we get 



(4.5) 

Taking derivative with respect to x in the both sides of the above relation leads to 

( e -*i(5)(x-6*) > if £ > ^ 

V ' (;r ' r) = j e^-)(l|a$ 1 ( ! ){ , - , ^( Z )e^) ) if < x < b*. (4 ' 6) 
It follows that V'(x, b*) < 1 when x > b*. Further, for < x < b*, 

V"(x,b*) = -Q^q^bW-*) (l + a^iq)^ V (?) (z)e- $l( '» z rfzj 

- a^^W^ib* -x) < 0. (4.7) 

Thus for < x < b*, V'(x, b*) > V'(b*, b*) = 1. 

Remark 4.1. From (4-6) and (4-7) we find that V'(x,b*) is continuous on (0, oo), 
V"(b*-,b*) = - a$i(g)W(«)(0) and \/"(6*+,6*) = -$i(g). So £/ia* if <J = 0, 

then V"(b*-,b*) ^ V"(b*+,b*) and, ifa^O, then V"(b*-,b*) = V"(b*+,b*). 

In summary, we have the following: 

Theorem 4.1. Consider the stochastic control problem (2.5). We have V(x) = V(x,b*) 
as defined in (4-5) and the threshold strategy with threshold b* is the optimal dividend 
strategy over S. 
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Abstract 

In this paper we study the optimal dividend problem for a company whose surplus 
process evolves as a spectrally positive Levy process before dividends are deducted. This 
model including the dual model of the classical risk model and the dual model with dif- 
fusion as special cases. We assume that dividends are paid to the shareholders according 
to admissible strategy whose dividend rate is bounded by a constant. The objective is 
to find a dividend policy so as to maximize the expected discounted value of dividends 
which are paid to the shareholders until the company is ruined. We show that the optimal 
dividend strategy is formed by a threshold strategy. 
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1 INTRODUCTION 



Recently, dividend optimization problems for financial and insurance corporations have 
attracted extensive attention. How should corporation pay dividends to its shareholders? 
A possible goal is that the company tries to maximize the expectation of the discounted 
dividends until possible ruin of the company. In recent years, quite a few interesting papers 
deal with the optimal dividend problem in the dual model. For example, Avanzi et al. 
(2007) consider the model which is dual to the classical risk mode, the authors studied 
how the expectation of the discounted dividends until ruin can be calculated when gain 
distribution has a exponential distribution or mixtures of exponential distributions and 
shows how the exact value of the optimal dividend barrier can be determined, if the jump 
distribution is known. Avanzi and Gerber (2008) examined the same problem for the dual 
model perturbed by diffusion. Moreover, they pointed out that "the optimal dividend 
strategy in the dual model is a barrier strategy. A direct proof that optimal strategy is a 
barrier strategy is of some interest but has not been given to our knowledge; the proof in 
Bayraktar and Egami (2008) is for exponential gains only." Yao et al. (2010) considered 
the optimal problem with dividend payments and issuance of equity in a dual risk model 
without a diffusion, assuming proportional transaction costs, they found optimal strategy 
which maximizes the expected present value of the dividends payout minus the discounted 
costs of issuing new equity before ruin. In addition, for exponentially distributed jump 
sizes, closed form solutions are obtained. Dai, et al. (2010, 2011) considered the same 
problem as in Yao et al. (2010) for a dual risk model with a diffusion with bounded gains 
and exponential gains, respectively. Avanzi, Shen and Wong (2011) determined an explicit 
form for the value function in the dual model with diffusion when the gains distribution is a 
mixture of exponentials. They showed that a barrier dividend strategy is also optimal and 
conjectured that the optimal dividend strategy in the dual model with diffusion should 
be the barrier strategy, regardless of the gains distribution. Until recently, Bayraktar, 
Kyprianou and Yamazaki (2013) using the fluctuation theory of spectrally positive Levy 
processes, show the optimality of barrier strategies for all such Levy processes. 
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Barrier strategies often serve as candidates for the optimal strategy when the dividend 
rate is unrestricted. However, if a barrier strategy is applied, ultimate ruin of the company 
is certain. In many circumstances this is not desirable. This consideration leads us 
to impose restriction on the dividend stream. Ng (2009) considered the dual of the 
compound Poisson model under a threshold dividend strategy and derived a set of two 
integro-differential equations satisfied by the expected total discounted dividends until 
ruin and showed how the equations can be solved by using only one of the two integro- 
differential equations. The cases where profits follow an exponential or a mixture of 
exponential distributions are then solved and the discussion for the case of a general 
profit distribution follows by the use of Laplace transforms. He illustrated how the optimal 
threshold level that maximizes the expected total discounted dividends until ruin can be 
obtained. In this paper we provide a uniform mathematical framework to analyze the 
optimal control problem with dividends payout for a general spectrally positive Levy 
process. We assume that dividends are paid to the shareholders according to admissible 
strategies whose dividend rate is bounded by a constant. Under this additional constraint, 
we show that the optimal dividend strategy is formed by a threshold strategy. This 
problem has been considered by Asmussen and Taksar (1997), Jeanblanc-Picque and 
Shiryaev (1995) and H0jgaard, B. and Taksar, M. (1999) in the diffusive case. 

The rest of the paper is organized as follows. Section 2 presents the model and for- 
mulates the dividend optimization problem. Section 3 discusses the threshold strategies. 
Explicit expressions for the expected discounted value of dividend payments are obtained, 
and in Section 4 we give the main results, it is shown that the optimal dividend strat- 
egy is formed by a threshold strategy. This strategy also called the refraction strategy, 
prescribes paying no dividends when the net surplus of the company is below an optimal 
barrier b*, and paying dividends at the fixed maximal rate a when net surplus exceeds b*. 
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2 The model and the optimization problem 



Let X = {X t }t>o be a spectrally positive Levy process with non- monotone paths on a 
filtered probability space (Q, J 7 , F, P), where F = (J r t)t>o is generated by the process X 
and satisfies the usual conditions. The Levy triplet of X is given by (c, a, IT), where c > 0, 
a > and II is a measure on (0, oo) satisfying 

POO 

/ (1 A x 2 )H(dx) < oo. 
Jo 

Denote by P x for the law of X when X = x. Let E x be the expectation associated with 

P x . For short, we write P and E when X = 0. The Laplace exponent of X is given by 

1 If 00 
^(6) = - \ogEe~ eXt = c6 + -a 2 6 2 + J (e' 9x - 1 + 9xl {lxl<1} )U{dx), (2.1) 

where 1a is the indicator function of a set A. In the sequel, we assume that —^'(0+) = 
E(Xl) > which implies the process X drifts to +oo. It is well known that if J" 00 yYi(dy) < 
oo, then E(Xl) < oo, and E(Xl) = — c + yll(dy). Note that X has paths of bounded 
variation if and only if 

POO 

a = and / (1 A x)U(dx) < oo. 

In this case, the Laplace exponent (2.1) can be written as 

1 f°° 

tf(0) = co + -a 2 9 2 + (e- 6x - l)U(dx), (2.2) 
^ Jo 

with c = c + J () 1 a;n(rfa;) the so-called drift of X. 

For an arbitrary spectrally positive Levy process, the Laplace exponent \P is strictly 
convex and limg^oo ^(9) = oo. Thus there exists a function $ : [0, oo) — > [0, oo) defined 
by $(g) = sup{# > O\V(0) = q] such that = g, g > 0. 

For more details on spectrally positive Levy processes, the reader is referred to Bertoin 
(1996) and Kyprianou (2006). 

Assume the canonical decomposition of X is given by 

X t = -ct + aB t + J t ,t> 0, (2.3) 



where {B t ,t > 0} is a standard Wiener process, {Jt,t > 0} is a pure upward jump 
Levy process that is independent of {B t ,t > 0}. In addition J — 0. Note that the 
dual model with diffusion in Avanzi and Gerber (2008) corresponds to the case in which 
U(dx) = \F(dx), where A > is the Poisson parameter and F is the distribution of 
individual gains, and the rate of expenses is given by Cq = c + xll(dx). In particular, 
when <t = 0, the model reduces to the so-called dual model in Avanzi, Gerber and Shiu 
(2007). 

We now recall the definition of the q— scale function W^ q \ For each q > there exits a 
continuous and increasing function : R — > [0, oo), called the g-scale function defined 
in such a way that W^(x) = for all x < and on [0, oo) its Laplace transform is given 
by 

poo 1 

/ e- ex W^(x)dx = — , 9 > *(q). (2.4) 

Jo - <i 

Closely related to W {q) is the function given by 

Z^(x) = 1 + q f W {q \y)dy, xER. 
Jo 

We will frequently use the following function 

Z (9 \x) = [ Z {q \z)dz, x G R. 



Note that 

Z (q \x) = 1, Z iq \x) =x, x<0. 

The following facts about the scale functions are taken from Chan, Kyprianou and Savov 
(2011). If X has paths of bounded variation then, for all q > 0, W^|(o,t») £ C 1 (0,cx)) 
if and only if II has no atoms. In the case that X has paths of unbounded variation, it 
is known that, for all q > 0, W^ (<?) |( ,oo) e C^O, oo). Moreover if a > then C^oo) 
may be replaced by C 2 (0, oo). Further, if the Levy measure has a density, then the scale 
functions are always differentiable. It is well known that 



W (5) {0+) = 



— , if X has paths of bounded variation, 



CO 



0, otherwise, 
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and 

[ if o- ^ 0, 

^( 5 )'(0+) = J i+moo) ^ if x ig compound Poisson 

1 c 

! oo, if a = and 11(0, oo) = oo. 

In all cases, if£pd) > 0, then W(oo) = l/E{X 1 ). If g > 0, then ~ e*^/*'^?)) 

as x — )• oo. 

We assume that the surplus process of the company is modeled by (2.3) if no dividends 
are paid. An admissible (dividend) strategy n = {L^\t > 0} is given by a nondecreasing, 
right-continuous and F-adapted process starting at 0. Let U n = {U^ : t > 0} be the 
company's surplus, net of dividend payments, at time t. Thus, 

U? = X(t)-L?, t>0. 

In this article we are interested in the case that ir only admits absolutely continuous 
strategies such that 

dDf = F(t)dt, (2.5) 

and for t > 0, l*(t) satisfies 

< F{t) < a, (2.6) 

where a is a ceiling rate. Moreover, we make the assumption that, if X has paths of 
bounded variation, 

a<c+ [ yU(dy). (2.7) 
Jo 

We define the dividend-value function V w by 



V n (x) = E 



l \- qt F {s)ds\Ul = x 
Jo 



where q > is an interest force for the calculation of the present value and tv is the time 
of ruin which is defined by 

r n = M{t > 0\U? = 0}. 

We denote by S the set of all the admissible dividend strategies. The objective is to solve 
the following stochastic control problem: the maximal dividend-value function, which is 
defined as 

V(x) = supV n (x), (2.8) 



and to find an optimal policy ir* G S that satisfies V(x) = V n *(x) for all x > 0. In this 
paper, we will prove that the optimal dividend strategy is formed by a threshold strategy 
with parameters b* (the definition of b* is given by (4.4)) and a: whenever the controlled 
risk process is below b*, no dividends are paid; however, when the controlled risk process 
is above this level, dividends are paid continuously at the maximal admissible rate a. 

3 Threshold dividend strategies 

In this section, we assume that the company pays dividends according to the following 
threshold strategy governed by parameters b > and a > 0. Whenever the modified 
surplus is below the threshold level b, no dividends are paid. However, when the surplus 
is above this threshold level, dividends are paid continuously at a constant rate a that does 
not exceed the rate of expense c. We define the modified risk process Ub = {Ub(t) : t > 0} 
by Ub(t) = X t — Db(t), where Db(t) = a J* * lsjj b {t)>b}dt. Let Db denote the present value 
of all dividends until time of ruin T, 



where T = inf {t > : U b (t) = 0} with T = oo if U b (t) > for all t > 0. Here q > is the 
discount factor. Denote by V(x,b) the expected discounted value of dividend payments, 
that is, 



Clearly, < V(x,b) < J and linr^ V(x, b) = |. 

Denote by T the extended generator of the process X, which acts on C 2 function g 
defined by 



^ Jo 
Theorem 3.1. Assume that, as a function of x, V(x, b) is bounded and twice continuously 
differentiable on (0, b) U {b, oo) with a bounded first derivative. Then V(x, b) satisfies the 




V(x,b) = E(D b \U b (0)=x). 
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following integro- differential equations 

AV(x, b) = qV(x, b), 0<x<b, (3.2) 

AV(x, b) - aV'(x, b) = qV(x, b) - a, x > b, (3.3) 
and initial condition V(0,b) = 0. 

Proof Applying Ito's formula for semimartingales one has 
E x [e-^V(U b (tAT),b)] = V(x,b) 

rthT 

+E X / e-«'[(A-q)V(U b (s),b) - al {Ub{s)>b} V'(U b (s),b)]ds. 
Jo 

Letting t — > oo and note that 1^(0, b) = we have 

V(x, b) = aE x / e- qt l {Ub[t)>b} dt 
Jo 

if and only if 

(A - q)V(x, b) - al {x>b} V'(x, b) = -al {x>b} . 

This ends the proof. 

Remark 3.1. From (3.2) and (3.3) one can prove that V satisfying the following 

V(b-,b) = V(b+,b)=V(b,b), (3.4) 
(ci + a)V'(b+, b) - dV'ib-, b) = a, if a = 0, (3.5) 
V'(b+, b) = V'(b-, b), if a ^ 0. (3.6) 

Define the first passage times, with the convention inf = oo, 

T+ = inf{t > : U b (t) > b}, T b = M{t > : U b (t) < b}. 

For q > 0, let 

^(q) = SU p{6 > 0\V(0) + a6 = q}. 

The following result generalized the result of Ng (2009, Theorem 2) in which only the 
dual of the classical insurance risk model was considered. 



Theorem 3.2. For x > b, 



V(x,b) = \ V(b,b) - ^ ) e-* 1 ^*-*). (3.7) 



Proof By using the strong Markov property of XJ h at T b as in Yin et al. (2013), we 
have 

V(x,b) = ^ - ^E x (e-« T >~ ,T~ < oo) 

+E x (e-^V(U b (T-),b),T- < oo). 

The result (3.7) follows since P x (Ub(T^) = b) = 1 and 

E x (e-^~,T- < oo) = exp(-$ 1 ( ? )(x - &)). 

Theorem 3.3. For < x < b, 

V(x,b) = --a 2 W {q \b-x) 
2 

+ V(b, b) ^feW-*) - i<7 2 $i(g)W^(6 - x) 

+a<$> 1 {q)e^ q){b - x) J* * W i - q \z)e~^ q)z dz^j 
+ ^ Q<r 2 $i(g)W (9) (& - ^) - e* l(9)(6 - x) $i( ? ) * Z (9) (^)e-* l(9)z ^ 

^ e *i(5)(b-x) ( ? _ a$1 ( g )) ^ X W^\z)e-^ q)z dz. (3.8) 



— t 

Q 



where 



Here 



V(b,b) = 



\<j 2 W^{b)- B(b)-C(b) 



A(b) 



A(b) = e $1 < 9 »C - la 2 ^ 1 (q)W^(b) + a$ 1 (g)e* 1 ^^ f W^(z)e-^ z dz, 

5(6) = ^ Q(7 2 $i( g )W (9) (&) - e* l( « )(6) $i(g) J" Z (q \z)e-^ q)z dz^ , 

C(b) = (g - a$i(g)) / W^e-* 1 * 9 ^. 

9 Jo 



In particular, when a = 0, 



V(x,b) = V(b,b)e* l{q){b - x) (l + a$i(g) / W {q \z)e-^ q)z dz 



b—x 



+ -e* 1 ^ 6 ^ (q - a$i(g)) " W (q \z)e~^ z dz 
Q Jo 



« e *i(9)(6-x) 



(3.9) 



where 



V(b,b) 



*Qi(q) Z^(z)e~^ q >dz - f (g - /^(^(^e"* 1 ^^ 
1 + | 6 W^(z)e-^(i>dz 



Proof We first assume that A := J °° Tl(dy) < oo. Set c = c + yll(dy) and 



f{y)dy 



, then / is a probability density on (0, oo). In this case, in view of (3.7), 



the integro-differential equation (3.2) can be written as 



l -a 2 V"(x, b) - coV (x, b) 



-X 



a 



V(b,b)-- 



exp(-$i(q)(x + y- b))f(y)dy 



b—x 



-A 



pb—x 

/ V{ 

JO 



Xa 



x + y,b)f(y)dy-—(l-F(b-x)) 



+ (A + q)V(x,b), < x < b, 



(3.10) 



where F is the distribution function of /. Replace the variable x by z = b — x, and define 
W by W(z, b) = V{b - z, b), < z < b. The (3.10) becomes 



-a 2 W"(z,b) + c W'(z,b) 



-A 



a 



W(0,b) - - 
Q 



/oo 
exp(-$i(g)(j/ - z))f(y)dy 



Xa 



-X W(y,b))f(z-y)dy (1 - F(z)) 

Jo Q 

+(X + q)W(z,b), 0<z<b, 



(3.11) 



with initial condition W(0,b) = V(b, b) and boundary condition W(b,b) = 0. We extend 
the definition of W by (3.11) to z > and denote the resulting function by w. Then we 
have 

1 



a 2 w"(z) + c w'(z) 



-X 



w(Q) - 



exp(-$i(g)(y - z))f{y)dy 

-X [ Z w ( y ))f( z - y )dy-—(l-F(z)) 
Jo 1 

+ (X + q)w(z), z>0, (3.12) 
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with w(0) = V(b, b) and w(b) = 0. 



For a function g, denoted by g the Laplace transform of g, i.e. g(£) = J °° e ^ y g(y)dy. 
Then the Laplace transform w for w can be easily determined from Eq. (3.11) as 

. M" 1 + fr(0)) + cow(0) + %(f(0 - 1) - F^y^W - f )(M(?)) - /(0 

= ; • 

\° 2 e + <tf - (A + q) + A/(0 

(3.13) 



Note that 



POO 1 

/ e _a *W (9) (a;Wa; = , — , 

A |^e + c e-(A + g) + A/(O 

/•oo /*a; i 

/ e'^dx W {q \y)dy = s , 

■/o 7o £(§<r^ + ^ - (A + q) + A/(0) 

/•oo fx py 

/ e-^dx / dy / W (<?) (z)d 
Jo Jo Jo 



1 



e(^ 2 e + c e-(A + g) + A/(o) 

Now inverting (3.13) gives 

«;(*) = - l -a 2 W^\z)+ l -a 2 w^) fw^\z-y)5' Q {y)dy 

+c o w(0)W^(z) + Xw(0) (W® * (f - f($i(q))5 ) * l) (z) 



\a ( Z iq) - 1 \ 
+y (*i(q)—j— *(6 -f)*l) (z) 

+^m i (q))-l)(W^*l)(z), (3.14) 

where 5o is the Dirac delta function at 0, hi * hi stands for convolution of hi and hi and 
l(z) = exp(&i(q)z). After some tedious calculations, we get 



W^(z-y)5' (y)dy = W^\z), 

X(w^*(f-f(^ 1 (q))6 )*ij(z) = a$i(q) J* W^\z - y)e^ y dy 

+e *i(5)* _ l a 2 W {q) \z) 
-(^i(q) + Co)W^(z), 

\(Z®-l)*(6o-f)*l)(z) = Q<7 2 g$i(g) + cog^ J* W^(z - y)e*^ y dy 

+ \a 2 qW {q \z) -q f Z {q \z - y)e^ q)y dy. 
2 Jo 
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Substituting the three expressions above into (3.14) we arrive at 



w(z) = -\a 2 W {q \z) 

+ w(0) (e* 1 ®' - ^a 2 ^ 1 (q)W^(z) 

+a$ 1 (g)e* l(<?)z jf W (q \y)e~^ y dy 
a ( 1 



+ ^ \^a 2 ^ 1 (q)W^z-e^ z ^ 1 (q) j (y^^dy 

+ - e *^ z (q-aQ^q)) f W {q \y)e~^ {q)y dy, 
1 Jo 

and the result (3.8) follows since V(x, b) = w(b — x) and w(0) = V(b, b). 



(3.15) 



Now, we assume that A := / °° H(dy) = oo. Let IT n be measures on (0, oo) defined by 

Tl n (dx) = n(dx)lr ( i i0on (x), n > 1. 



Then we have 



U n (dx) < n 2 J x 2 Yi(dx) + J (1 A a; 2 )IT(c£r) 



< oo. 



Set c n = c + Jq yll n (dy) and f n (y)dy = then for each n > 1, / n is a probability 

density on (0, oo). Similar to (3.10) we consider the following integro-differential equation 

1 



^a 2 V n "(x,b)-c n V n '(x,b) = -X 



a 



V n (b,b)-- 
1 



exp(-9 n (g)(:r + y- b))f n {y)dy 



b—x 



rb ~ x Xcv 
-\ I V n (x + y, b)f n {y)dy (1 - F n (b - x)) 



io q 

+(X + q)V n (x,b), < x < b, 



(3.16) 



where Q n (q) = sup{6> > 0\^ n (9) + ad = q}. Here, 

1 f°° 

*n(0) = cj + -a 2 6 2 + J (e' 9x - 1 + 9xl {lxl<1} )U n (dx). (3.17) 
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Repeating the same argument as the case that A := f£°H(dy) < oo, we obtain 
V n (x,b) = - l -a 2 W^\b-x) 

+ V n (b, b) ( e @ ^ b ~^ - l -a 2 Q n (q)W^(b - x) 



2 

"b—x 



+ae n (q)e @n{q){b ~ x) / W { a q) (z) e - @n{q)z dz 
Jo / 

+ ^ (^a 2 Q n (q)Wi q \b -x)- e &n(q)ib - x) e n (q) J* * Z {q) (z)e- 0n{q)z dz^j 

+ -e 0n{q){b ~ x) (q-ae n (q)) f X Wi q \z)e- 0n{q)z dz, (3.18) 
Q Jo 

where W^f 1 and Z$ are scale functions corresponding to the process X n with Laplace 
exponent \I/ n . Since lim^oo ^ n (9) = ^(9), then X n — > X weakly as n — > oo. Thus 
lim^K = V, lim^W^ = Wb\ lim^ = and lim^ Q n = $ x . Conse- 
quently, (3.8) still holds for this case. This ends the proof of Theorem 3.3. 

Remark 3.2. For a > 0, from the graph of 

#($i(g)) + a$i(g) = q 

one can verify that $i(g) — > when a — > oo. After some simple calculations we get 
]im a -> 00 a$i(q) = q, 

lim a(q - a®i(q)) = qV'(0+), 

and 

V T 7-/1 ,V ^(fc) ^(0 + ) *'(0 + ) 

hm V(b, b) = — T77T + -7 , ./■ -. 

«^oo v ' ; Z(«)(6) gZ(«)(6) g 

As a result, for < x < b, we arrive at 

J5t^)-£|^-.)-^-.) + ^(^-i). 

which is the expected discounted value of dividend payments for the barrier strategy. See 
Lemma 2.1 of Bayraktar, Kyprianou and Yamazaki (2013). 

4 Optimal dividend strategy 



Suppose the maximal dividend- value function V is C 2 (0, oo) (resp. C 1 (0, oo)) when X is 
of unbounded (resp. bounded) variation. Standard Markovian arguments, see Fleming 
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and Soner(1993), formally yield that V(x) satisfies the following Hamilton- Jacobi-Bellman 
(HJB) equation 

TV(x) - qV(x) + sup {r(l - V'(x))} = (4.1) 

0<r<a 

with V(0) = 0, where T is the extended generator of the process X, which defined by 

1 f°° 
Tg(x) = -a 2 g"(x) - cg'(x) + J [g(x + y) - g(x) - g'(x)yl {lyl<1} ]U(dy). (4.2) 

It follows from the HJB equation (4.1) that an optimal dividend policy has to fulfil r = 
if V'(Uf_) > 1; r = a if V'(U^_) < 1. In some situations the optimal dividend strategy 
is a threshold strategy. In fact if V'(x, 0) < 1 for x > 0, then the threshold strategy 
with b* = is optimal, if V'(x, b*) > 1 for x < b* and V'(x, b*) < 1 for x > b*, then the 
threshold strategy with b* > is optimal. The optimal threshold b* can be obtained by 
V'(b* : b*) = 1. From those facts one sees that if V(x, b*) is a concave function on (0, oo), 
then the optimal dividend strategy is a threshold strategy. 

It follows from (3.7) that 

V'(x, b) = (v(b, b) - ^ e -*i(«X*- b ), x > b. (4.3) 

If $i(g)f < 1, then V'(x,0) < 1 since F(0,0) = 0. Thus b* = 0. Now, suppose that 
$i(g)^ > 1, we get the condition for b*: 

-$!(?) ^(6*,6*)-^=l. 
Or, equivalently, 6* is the solution of the equation 



From (3.7) and (3.8) we get 



(4.5) 

Taking derivative with respect to x in the both sides of the above relation leads to 

, e -*i(«)(*-&*), ifx>6*, 
)-< e *,( 5 )(i'-x) ^i + a^g) f-^M^e-*'^^) , if0<x<b*. (4 ' 6) 
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It follows that V 



(x, b*) < 1 when x > b*. Further, for < x < b* , 



V"(x, b*) 



-$ 1 (g)e <E>l(9)(b *- :c) 1 + / W {q \z)e 




^ z dz 



) 



a$i(g)W (9) (6* -x) < 0. 



(4.7) 



Thus for < x < b*, V'(x, b*) > V'(b*, b*) = 1. 

Remark 4.1. From (4-6) and (4-7) we find that V'(x,b*) is continuous on (0, oo), 
V"(b*-,b*) = - a$i(g)W" (9) (0) and V"(b*+,b*) = So that if a = 0, 

then V"(b*-,b*) j£ V"(b*+,b*) and, ifa^O, then V"(b*-,b*) = V"(b*+,b*). 

In summary, we have the following: 

Theorem 4.1. For any spectrally positive Levy process (2.3), consider the stochastic 
control problem (2.8). Let S be the class of admissible dividend strategies satisfying (2.5)- 
(2.7). Then we have V(x) = V(x,b*) as defined in (4-5) and the threshold strategy with 
threshold b* is the optimal dividend strategy over S. 

Remark 4.2. It is interesting to note that the optimal strategy with bounded rate of 
dividend payment is formed by a threshold strategy for spectrally positive Levy process 
regardless of the Levy measure. However, for the spectrally negative Levy model —X, it 
was shown in Kyprianou, Loeffen and Perez (2012) that the optimal strategy is formed by 
a threshold strategy under certain regularity condition on the Levy measure of —X . 
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